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1 Introduction 

The question about the existence of absolutely stable strange quark matter in astro- 
physical compact objects is still debated in a number of theoretical and observational 
investigations, see Ref. [1] for a recent review. In particular the formation of quark 
stars is a very interesting issue. In the seminal paper [2], the process of conversion of 
a neutron star into a quark star was shown to be a very exothermic process which, 
if it really occurs in Nature, should lead to powerful explosive phenomena, maybe 
also connected with supernovae and gamma-ray-bursts lasiE]. The relevance of this 
process, from the phenomenological point of view, is related to the time scale of the 
conversion. Only if sufficiently fast (with time scales of the order of tens of seconds 
or less) the formation of a quark star could release detectable signals; on the other 
hand, a very slow conversion would not probably provide any evident signature. In¬ 
terestingly, semi-analytical estimates and numerical hydrodynamics simulations have 
provided hints for a very fast conversion, occurring on time scales of the order of 
ms pill El El [To]. In this short contribution, we will review the approach used for 
studying the conversion of a neutron star into a quark star based on the assumption 
of a inhnitely thin combustion zone and we will discuss why, in this scheme, the 
combustion stops before the whole hadronic star is converted. 

2 Combustion of hadronic stars 

The process of combustion is a very complicated phenomenon which in principle 
must be modeled by coupling the hydrodynamic equations for the mixture of the 
two fluids (the fuel and the ashes) with the equations of conservation of chemical 
species in which both diffusion processes and rates of chemical reactions are included 
HU, see [121 E] for the case of quark stars. A common assumption, especially used 
in astrophysical problems such as type la supernovae [T3|, consists in treating the 
combustion zone as a surface of discontinuity separating the fuel from the ashes. 
This assumption is necessary since in many cases the width of the combustion zone is 
much smaller then the size of the system and it would be computationally unfeasible 
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to resolve the microscopic dynamics of the combustion zone. As in the case of shock 
waves, in presence of a discontinuity, hydrodynamic equations lead to the conditions 
of continuity of the baryon flux, the momentum flux and the energy flux across the 
burning front [g. Indicating with Cj, pi, rii, Wi, Xi = (e^ +Pi)l'nl the energy density, 
pressure, baryon density, enthalpy and dynamical volume of the i-th fluid, those 
conservation laws, generalized to relativistic hydrodynamics |T 5 ], read: 


riiUi 

= n2U2 = j 

(1) 

(P2-Pi)/(Ai-X2) 

= f 

(2) 

m;i(Pi,Xi)Xi - W2 {p2 ,X2)X2 

= (Pi-P2)(Ai+X2) 

( 3 ) 


where Ui are the four-velocities of the two fluids in the front rest frame and j is the 
baryon number flux across the front and it must be determined within a microscopic 
kinetic approach, as the one of Ref. [ 2 ]. This system of equations allows to determine 
the velocities of the two fluids and the state of fluid 2 once the initial state of fluid 1 
is hxed. Depending on the values of the velocities Ui and the sound velocities of the 
fluids Cj, one can obtain detonations, which are processes of combustion driven by 
a shock wave, or deflagrations in which instead combustion proceeds thanks to the 
diffusion of heat or of chemical species m. 

The condition of exothermic combustion, generalized to relativistic hydrodynamics 
by Coll in Ref. [ 15 ], reads A(p, X) = ei(p, X) — e2(p, X) = Wi(p, X) — W2{p, X) > 0 
and allows to hnd the window of baryon densities of fluid 1 for which the combustion 
can proceed. We want here to clarify the meaning of this condition. Let us £x the 
initial state A of fluid 1 : pi = PaXi = = Xa (the temperature of fluid 1 

is set to zero). From this point one can draw the shock adiabat of fluid 1 in the 
(p,X) plane. Moreover, by using the equation of state of fluid 2 and Eq .3 one can 
draw in the same plane the detonation adiabat, see Fig.l. In the following we will 
demonstrate that, the Coil’s condition implies that the shock adiabat lies below the 
detonation adiabat. Let us assume that the equation of state of fluid 2 is a generic 
polytrope 62 = cm2 +^2/(7 — 1) ([E]), P2 = knl, where 7 is the adiabatic index and 
1 < 7 < 2 (the second inequality implying that the equation of state is causal at all 
densities ini)- 

One can derive the following expression for the energy density as a function of p2 
and X2: 


62(^2, ATa) 


( 3:^(7 - 1 ) + 2P2X2 + 

2X2(7 - 1) 


( 4 ) 


Let us £x Xi = X2 = Xa and assume A{pa,Xa) > 0 . If P2 > Pi = Pa, fhe 
detonation adiabat lies above the shock adiabat. With this setting, the detonation 
adiabat reads (adding and subtracting e2{pA, Xa))- 


A{pa, Xa) =PA-P 2 + 62 (^ 2 , Xa) - e 2 {pA, Xa) 


( 5 ) 
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Figure 1: Illustrative plot of the shock and the detonation adiabats in the case in 
which the Coil’s condition is fulhlled or not. 

which after some manipulation and using Eq.4 reads: 

A{pa,Xa) = - 1) + 47P2^a - ^^^(7 - 1) + (6) 

+ (P 2 -Pa)^ (7) 

7-1 

Since 1 < 7 < 2 the sign of A(p^, Xa) clearly determines the sign of P 2 —PA- Thus, if 
A{pa, Xa) > 0, i.e. if the Coil’s condition holds true, the initial point A lies below the 
detonation adiabat 0. Since the detonation adiabats and the shock adiabats do not 
cross in the (p,X) plane [H] this implies that the whole shock adiabat lies below the 
detonation adiabat (at least in the standard cases).The Coil’s condition is necessary 
for obtaining detonations. Detonation is a process of combustion which is driven by 
a shock wave propagating within the fuel: when the shock wave passes through A, 
the fluid is compressed and heated up to the state A* which lies on the shock adiabat 
(see Fig.l). The chemical reactions start, the fluid expands and cools down until the 
combustion is complete and the state B, lying on the detonation adiabat, is reached 
[13]. On the other hand, if the shock adiabat of the fuel lies above the detonation 
adiabat, it is not possible to trigger the combustion via a shock wave. A shock wave 
passing through the point A would heat up the matter but not at a sufficiently high 
temperature to start the combustion (see Fig.l). Once the Coil’s condition is fulhlled 
the hydrodynamical combustion can proceed either as a detonation or a dehagration 

^For values of 7 > 2 , one can still obtain the same conclusion but depending on the specific values 
of a, PA and Xa 
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depending on the specific microphysics and on the equation of state. The equation 
A{pa,Xa) = 0 allows to find the state A, characterized by its baryon density 
below which the combustion cannot proceed anymore. At this density ei = 62 = e^, 
Pi = P2 = Pa which together with Xi = X2 = Xa also implies ni = n2 = n\. At 
this value of baryon density there is no surface of discontinuity anymore and the two 
phases are in mechanical equilibrium. Previous studies have shown that, for many 
equations of state, n\ ~ 0.2 — 0.3 fm“^ [6l [ 9 ] and indeed the fast combustions found 
in the numerical simulations of Ref. | 9 ] stop exactly at those values of density. 

One has to remind that the scheme based on a infinitely thin combustion zone is 
only an approximation. Relieving this approximation and considering the microscopic 
processes of diffusion of quarks and reactions between quarks (such as u + d— > u + s) 
occurring in the finite width combustion zone allows to follow the subsequent evolution 
of the system which proceeds until the whole star is converted [T8] . 

3 Conclusions 

We have discussed the approximation scheme of combustion based on the assumption 
of a infinitely thin combustion layer. This model has been widely used in numerical 
simulations of type la supernovae [TH] and recently also for numerical investigations of 
the conversion of hadronic stars into quark stars At densities larger than the 

critical density n\, for which the Coil’s condition is fulfilled, the process of conversion 
proceeds very fast with the effective velocity of conversion significantly augmented by 
Rayleigh-Taylor instabilities. After a few ms, when the conversion front reaches n\, 
a big part of the star is converted but a few O.IMq remain unburnt and will convert 
on a longer time scale, of the order of tens of seconds [IH]. Note that, up to now, 
the conversion process has been studied only within cold and non rotating hadronic 
stars. The dynamics of the birth of quark stars can be qualitatively different in other 
astrophysical situations such as supernovae and protoneutron stars and mergers of 
neutron stars. Those cases are, to date, essentially unexplored and detailed investiga¬ 
tions are therefore needed to better clarify the scenario of coexistence of two families 
of compact stars, hadronic stars and quark stars, proposed in [ 191 1 ^ - 
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